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SCHRODINGER EQUATION FOR THE CHIRAL 
GEOMETRY IN ATOMIC NUCLEI 


Radu BUDACA! 


Abstract. The dynamical and spectral features of the chiral geometry in atomic nuclei are 
described by means of a triaxial rigid rotor Hamiltonian cranked by quasiparticle 
alignments. The problem is treated alternatively in the space of angular momentum states 
and using a Schrédinger equation for a continuous variable associated with an angular 
momentum projection. The later is constructed using a semiclassical approach. Numerical 
applications performed for various deformation and alignment conditions with both 
methods are compared in terms of energy levels and wave functions. 
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1. Introduction 


The concept of chirality or handedness is related to the asymmetry properties of a 
system. Basically, a chiral system is distinct from its mirror image. This property is 
systematically encountered in nature and subsequently is important in most do- 
mains of science. One of the most immediate examples of chirality are the human 
hands. Shells of snails, spirals and coils of plant veins, or flower formations are 
other macroscopic biological ensembles with a chiral geometry. From mathematical 
point of view, a system is chiral if its mirror image cannot be obtained by a 
combination of only translation and rotation transformations. In science, the 
chirality is more often invoked in molecular and biological chemistry, where chiral 
partners of compound molecules often exhibit drastically different properties such 
as smell, taste, or colour. A notable example of a chiral biological molecule, is the 
DNA double helix. The chirality in physics is mostly associated with spin dynamics 
and the polarization of the electromagnetic waves. 


In nuclear physics, chirality is related to the trihedral geometry of three mutually 
perpendicular angular momentum vectors corresponding to the core rotation and 
two distinct spins coming from out of core nucleons. The three vectors can therefore 
be arranged into a right-handed or left-handed order. Given the fact that vectors are 
involved instead of coordinate loci, the nuclear chirality is of the dynamical type. 
The two chiral configurations in this case can be obtained from each other through a 
rotation and a time reversal transformation instead of a space inversion employed for 
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the static chirality. The conditions in which the three spin vectors achieve a mutual 
orthogonality were originally proposed by Frauendorf and Meng [1]. It all comes 
down to the quasiparticle nature of the nucleons which couple to the triaxial core. 
More precisely, a nucleon of particle (hole) nature will rotate around the short (long) 
axis of the mean field generated by the triaxial core, which in turn favors the rotation 
around its medium body-fixed axis with the maximal moment of inertia (MOI). The 
schematic representation of such a geometrical arrangement is given in Fig. 1. 


A3(m) 


Th 


Figure 1: Schematic representation of the chiral geometry in atomic nuclei. Right- (R) and left- 


handed (L) configurations of rotor (R) and total (I ) angular momentum vectors are correspondingly 
labeled. For the direction of the total angular momentum vector one also pointed out the 
stereographic variables g and 6. The considered departure from the perfect orthogonal alignment of 


the single particle spins (7) and Gi ) is also presented. 


The two chiral configurations can be related either by the interchange between the 
quasiparticle alignments or considering an opposite rotation for the triaxial core. 
The distinction between the two opposite chirality configurations is valid in the 
body fixed frame of reference, but are indistinguishable in the laboratory reference 
system where the measurements are taken. Due to the restoration of the chiral 
symmetry, a degenerated pair of AJ = 1 rotational bands of the same parity should 
be observed. The first observation of such bands and consequently the experimental 
realization of chiral geometry was reported in a few N = 75 odd-odd nuclei [2]. This 
event triggered a concerted theoretical and experimental effort for establishing the 
extent of this dynamical mode in other nuclei. As a result, chiral bands are presently 
identified in more than fifty nuclei [3, 4], be it in its originally proposed one proton 
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and one neutron configuration, or with multiple quasiparticle alignments. Most of 
the experimental realizations can be localized within the so called islands of 
chirality centered on the mass numbers 4 = 80, 100, 130 and 190 with specific 
single-particle states involved in the realization of the chiral geometry. The 
quasiparticle alignments are usually generated by protons and neutrons from 
intruder orbitals go, A112, and i)32, assuring thus a high spin contribution to the total 
angular momentum of the system. The topic of nuclear chirality evolved also to a 
point where it is considered in connection to other nuclear structure properties or 
phenomena such as octupole correlations [5], shape phase transitions [6] and shape 
coexistence phenomena [7], or softness of the triaxial deformation [8]. 

The expected degeneracy of the chiral bands is rarely experimentally realized, being 
confirmed only approximately for few nuclei in limited ranges of total angular 
momentum. Nevertheless, the chiral assignment of these bands persisted due to the 
same underlying structure. At low spins, the core rotation is slow and the two chiral 
configurations are attained as turning points for the so called chiral vibration [2]. 
The excited state of the same angular momentum and parity is shifted in this case 
with a vibrational energy. With the increase of core’s rotation, the chiral directions 
of the total angular momentum stabilize, marking thus the beginning of the static 
chirality regime. The lack of degeneracy between the partner bands in this high spin 
regime can be ascribed to the deviations from perfect orthogonality of the all three 
vectors involved in the chiral geometry. Indeed, the uniform rotation of the core is 
plagued by precession components coming from its triaxiality [9, 10, 11, 12, 13], 
while the fractional occupancies of the coupled quasiparticles induce tilting to their 
assumed axial alignments [14]. Of course, in this situation the system is no longer 
chiral, but valuable information can be extracted having this picture in mind. 
Although the obvious theoretical tool to approach the chiral geometry in atomic 
nuclei is the Particle-Rotor Model (PRM) [15], its interpretation regarding the 
direction of the total angular momentum in relation with the nuclear density 
distribution was achieved using the Tilted Axis Approach (TAC) [16]. On one hand, 
the quantum nature of the PRM offers spectral signatures and selection rules for the 
electromagnetic properties of the chiral partner bands [17, 18, 19, 20, 21, 22]. On 
the other hand, the semiclassical upbringing of the TAC formalism is able to unveil 
the dynamical behaviour of the chiral system [23, 24, 25, 26]. Advances in 
understanding and describing the properties of chiral bands were driven by various 
approximations and extensions of these two models such as TAC plus random 
phase approximation [27, 28], collective Hamiltonian [29, 30, 31], boson expansion 
[32, 33], projected mean-field [34, 35, 36], and semiclassical models [9, 10, 37, 38]. 
Alternative theoretical approaches to those based on PRM and TAC such as the 
Interacting Boson Model [39, 40, 41] must also be mentioned. 

The semiclassical approaches deserve a special attention due to their ability to 
combine the advantages of the PRM and TAC formalisms. This is achieved by 
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ascribing to the original quantum PRM Hamiltonian a classical picture in terms of 
directional angles of the total angular momentum vector, from which the dynamics 
of the system can be better understood. The semiclassical method can also be used 
to recast the original quantum PRM Hamiltonian into a Schrédinger equation for a 
suitably chosen continuous variable, which can calculate experimental observables 
and relate them with the underlying classical motion. The same program had great 
success in the description of the wobbling excitations [42, 43, 44, 45, 46], a 
phenomenon having many common features with the nuclear chirality. To keep the 
formalism tractable, a frozen approximation is usually adopted for the quasiparticle 
spins, which actually relates the PRM and TAC approaches. Indeed, with adoption 
of rigid or frozen quasiparticle alignments, the PRM Hamiltonian becomes a 
cranked rotor Hamiltonian with constraints defined by the quasiparticle alignments. 
Such a model is then well suited for the investigation of realistic arrangements of 
the implicated spin vectors deviating from the perfect chiral geometry. The 
analytical formalism behind this quantum model is presented in Section 2. The 
construction of the Schrédinger equation associated to this quantum problem, by 
means of a semiclassical procedure, is detailed in Section 3. Given the fact that in 
order to obtain a Schrédinger equation, one must adopt additional approximation 
and perform a requantization, it is interesting to compare the eigensystems of the 
two approaches. This is the main focus of the study and is realized in Section 4 with 
various numerical applications. The conclusions of this comparison are synthesized 
in the last section. 


2. Constrained triaxial rotor model for chiral geometry 


The interaction between the rotation of a triaxial core and two quasiparticle (7) and 
(7'), is usually described by an appropriate extension of the PRM: 


H= Ar+ Hyp + Hp (2.1) 
The first term 


Hp= >> Ar(le—ie—-J'x)” (2.2) 
k=1,2,3 

is the triaxial rotor Hamiltonian for the core angular momentum (R)= (I )- qi )- G’ ) 
, with (I ) being the total spin of the system. Each principal axis rotation component 
is factorized by the corresponding inertial parameter Ax; = 12%), which is related 
to the moment of inertia 


4 2 
Tk = 30 sin” (> = jkr), (2.3) 
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defined in the hydrodynamic model [15] with the help of the triaxiality measure . 
The shapes of the triaxial core repeat themselves at y = 2/3 intervals, with 
redistributed axis assignments. The correspondence between the MOI and the 
core’s shape can be tracked with the help of the semi-axis lengths: 


1+ \ ZB cos (, - *1)| (2.4) 


where / is the axial deformation. This allows the identification of the intrinsic axes 
as the long (/), short (s) and medium (m) depending on the value of the triaxiality 
measure y. A triaxial rigid body will favor rotations around the medium intrinsic 
principal axis whose MOI is maximal. For the purpose of the present study, the 
third axis is chosen as a quantization axis and associated with the rotation of the 


Ry = Ro 


core. The third axis become medium for y € (60°, 120°), implying Ro < R3 < Ri, 
such that the other two axes 1 and 2, become long and respectively short. The 
alignment of the single-particle spins is determined by their quasiparticle character 
[47] as follows: a hole quasiparticle will rotate around the / axis, a particle around 
the s axis, while a nucleon in the vicinity of the Fermi surface will prefer rotations 
around the m axis. 

In what follows, the quasiparticle spins are considered rigidly aligned in respect to 
the body-fixed system of reference. This condition, also referred to as the frozen 
alignment (FA) approximation, is realized by replacing the quasiparticle spin opera- 
tors with their expectation values associated with the assumed alignment geometry. 
The perfect chiral geometry emerges when one quasiparticle is fully aligned to the 
s axis and another one to the / axis, with the core favoring rotations around the m 
axis. To encompass a more realistic situation, the spins of the quasiparticles are 
considered rigidly aligned along the principal planes 1-3 and 2-3 as in Fig. 1, with 
a. and a angles accounting for the deviations from the full / and s axial alignments: 


jr = ji, cosa, j4 =), 
jo = 0, j'g =j' cosa’, (2.5) 
jg = jsina, j’s =j' sina’. 


The FA assumption reduces the relevant part of the PRM Hamiltonian to [37]: 


a 


Boi. = AI? + Agi? + A3i? — 2A1jh; cosa — 2Aj' I cos a’ 

—2A3(j sina + j’sina’)I[3, i 
where one omitted the constant quasiparticle contribution. This operator defines a 
cranked rotor Hamiltonian for a total spin [, whose constraining terms originate 
from the quasiparticle alignments. In this way, the single-particle degrees of 
freedom are downgraded to perturbing effects for the total angular momentum. The 
Hamiltonian A «ig, can be diagonalized for each angular momentum value J, in a 
basis of 27+ | rotation states |[MK), where M and K are projections in the laboratory 
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and respectively body-fixed reference frame of the f;. The matrix elements of the 
A aign Can be expressed as: 
(IMK | Hatign|IM K') = (2.7) 


a (1? +1 — K®) + A3K (K — 2jsina — 2;' sin a’)| one 


+a’ [Vi+K +2)0—K -NT+K +10 K)ox,«--2 

+ \U- K+ HEF K—NT- KF T+ K)bK, 442] 

—Ajj cosa [Vi+ K+ T= Kx, + VI=K = T+ K)5x,0741| 
—i Aaj’ cosa’ [VE+ KFT) ox,40--1 = VI=K-DNT#R)5 x74] . 


The corresponding wave function is 
I 
IS:us) = > > Arxs|IKM), (2.8) 
K=-I 
where s is the solution’s order and Aj;x; are the corresponding eigenvector 
components. 


3. Schrédinger equation for a projection variable 


In order to obtain an alternative quantum description of the dynamics associated 

with the A aie Hamiltonian, one must pass through a semiclassical description. This 

stage introduces a new restricted set of complex variables, which will be ultimately 

quantized. The semiclassical description is based on a time-dependent variational 
0 


principle [48] ; 
sf (v a wat = 0, (3.1) 


applied to the A azign Hamiltonian with an appropriately chosen variational state. The coherent state 
for the angular momentum operators is the standard choice for the variational state [9, 10]: 


I+a\! [Pe giv7_ 
Wirar(zs¢)) = ( * ) eV ee" rygry 


Hatign — 


= s Spl OR +2) - 2) ete Our K). 


(3.2) 


Here one used the stereographic parametrization of the complex coordinates in 
terms of the azimuth angle g and a polar projection variable x = J cos 0 of the total 
angular momentum vector on the third (quantization) axis. Solving the variational 
principle, one obtains a classical counterpart of the quantum Hamiltonian A aien, 
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represented by an energy function depending on the stereographic variables of the 
variational function: 


H(z, y) = (Yim (2, ~)|Hatign|¥ru (2, ~)) 
(21 — 1)(I? — 22) 


= Ti (A; cos? y + Ag sin? y — A3) (3.3) 


—2\V/ I? — x?(A1j cosacos y + Aoj’ cosa’ sin y) 
I 
—2A3z(j sina + j’ sina’) + 5 (At + Ag) + AsI?, 


where f = | is assumed. The other product of the variational principle is a pair of 
equations of motion, 


{H,x}=%, {H,p} =¢, (3.4) 


which are expressed with the help of the Poisson bracket. Note the canonical form 
of the equations of motion, which is due to the chosen variables. Consequently, 
these variables are canonically conjugated {g, x} = 1, and one can assign x as being 
the generalized momentum, while ¢ as the generalized coordinate. This property is 
further used to quantize the classical energy function by means of a correspondence 
principle between the canonical variables and their differential operator 
counterparts. The procedure is however not straightforward, because the classical 
energy function have products of mixed rational and trigonometric functions of 
momentum and coordinate. In order to pass this obstacle, one recounts the 
properties of the classical energy function. It can be thus shown that the classical 
energy function (3.3) always has a single minimum in g when y is restricted to the 


(60°,120°) interval. The classical energy function can be then approximated by 
means of a harmonic expansion around the corresponding minimum points in g(x) 
for fixed values of x: 


~ 1 (PH 2 - 
H(z, p) * H(z, po(z)) + 5 | aa [yp — yo(z)]”. (3.5) 
2 \ 99" 7 go(2) 

With this, the trigonometric functions are replaced by powers of g, up to the second 


order. The minimum path function go(x) is defined as the solution of a quartic 
equation for cos@o(x): 


(22 — 1)V I? — 2?( Ag — A) cos Yo(x) sin yo(r) 
= 2I [Aoj’ cos g(x) cosa’ — Ajj sin yo(x) cos a] , (3.6) 


which originates from the minimum condition OH(x,g)/0g = 0. The role of the 
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harmonic approximation is twofold. Indeed, the classical energy approximated as a 
second order function in g, can now be quantized into a Schrédinger equation in 
the momentum space. Symmetrizing the mixed products of coordinate and 


a .d ‘ 
momentum, and then applying the correspondence principle g = I one arrives at 


the following differential equation: 


1 dF(z) B'(z) iyo(x)\ dF (x) 
 2B(x) dx? (aoe aay) dx * 


[¥0(x)]? — igo (a) 4 iyo(x)B' (a) 
2B(z) 2[B(x)]* 


pt yo(x)) + 


B"(z) _ [B(a)? 
4[B(a))?  2[B(x)/ 


F(z) = EF(z), (3.7) 


For an unknown wave function F(x). Prime and double prime denote x derivatives. 
Making the change of function 


f(x) = [B(x)-4 F(a)etS eee, (3.8) 
one can see that 
OPH(z, ~) a 
ae) = [eGo (3.9) 
‘ Og? J go(2) 
= are — x”)(Ag — A;) cos 2yo(x) + 2V/ I? — x? 


x [Aig cos Yo(z) cos a + Aoj’ sin (yo(x) cos a’) — ; 


plays the role of a coordinate dependent effective mass for the Schrédinger 
equation: 


-}_1_£ _ 1 _4 vq} g(a) = Ef (a). (3.10) 


2 / B(x) dx V B(a) dx 


| 1 1d 1 4d 
The effective mass defines the potential 


B'(x) — 9[B'(a))? 


8[B(a)P 32 (Ba) oe 


V(x) = H(x, o(x)) + 


and is also involved in the normalization condition for the new wave-function 
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I 
i \f (x) |? /B(x)dx = 1. (3.12) 


The Hamiltonian of Eq.(3.10) is similar to those constructed based on TAC 
formalism [29, 30, 31], but is completely determined just by the deformation and 
the alignment geometry of the treated system. 

If go(x) is known, then the two quantities B(x) and V(x) can be analytically 
described. When y = 90¢ (A; = A2 = A1,2), Eq.(3.6) has a simple solution, represented 


by a constant tan Dp) = ae Moreover, the potential V(x) acquires in this case a 


relatively simple compact expression: 


1 x? 
= 4/12 = op? Acaa ——— 
V(x) =vuo + 12 + veV I? — 2? + v3n* + U4 — eg 372° (3.13) 
where 
I 2I-1 
vp = Ail? + 3 (As — Ai.2), “a= (A3 — A1,2) , 
v1 = —2A3 (jsina + j’sina’) , v4 = -2, (3.14) 
1 
v2 = —2A12V/j2 cos? a + j? cos? a’, v5 = a. 
The corresponding mass function is also simplified: 
1 
B(z) = (3.15) 


\/2(I? — x?) Aj.2(j cosa + j/ cosa’) 


Note that B(x) is invariant to the change x @ —x, while V(x) is not. This aspect is 
actually generally valid, regardless of the triaxial deformation y. 


However, for a y # 90°, Eq.(3.6) must be solved numerically. Only one out of its 
four solutions is real and physically relevant. A continuous and differentiable 
instance of @o(x) to be used in the above definitions of the effective mass and 
potential is obtained as follows. The selected solution is calculated for each angular 
momentum in an extended array of values for x. As can be seen from the Eq.(3.6), 
the minimum condition is symmetric to the change +x, therefore one can limit the 
calculation only to x € [0, J). These points are then interpolated with a spline 
function. Few instances of the go(x) are shown in Fig.2. It is clear that the effect of 
the tilting becomes important only at high angles, being materialized in an 
extension of the functions’s range. 
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Figure 2: Solution go(x) of Eq.(3.6) for 7 = 7’ = 11/2 and selected values of the 
triaxiality parameter y = 80°, 90°, 100°, angular momentum / = 10 and 16, with and 
without tilting. 


The natural boundary |x| < J of the problem [49] leads to the inevitable choice of 
particle in the box wave-functions: 


file) = F(a) (B(a)\~4 (3.16) 
_ [B(x)]~4 ae _[ (Qn —1)aa . . | 2ntr 
= BT FS Aton [SDE 4 SY ain | aI |}. 


or the diagonalization basis. These functions satisfy the Dirichlet boundary 
condition F7 (J) = F7(-J) = 0. 

The basis dimension is truncated when a satisfactory convergence of eigenvalues is 
achieved. Previous numerical applications of the model established that a hundred 
basis states assures a satisfactory convergence of the diagonalization results for all 
relevant deformation and alignment conditions. The order s of the diagonalization 
solutions distinguishes the energy states belonging to distinct chiral bands. For a 
meaningful comparison with the exact diagonalization results, one expresses the 


total wave function as an expansion in rotation states 
I 


|Yias) = D> Arxs|\ZKM). (3.17) 
K=-I 


The expansion coefficients 
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2 
2 


F#(K), (3.18) 


I 
Alks = | So pi(K) 


K=-I 
are determined with the help of the continuous probability distribution 
pi (a) = |Fi(a))” = |fi(@)? VBQ), (3.19) 
resulting from the diagonalization of the Schrédinger equation. 


4. Numerical application 


The energy splitting between the yrast and yrare partner bands 
AE(I) = Fs=9(I) — Es=1(1), (4.1) 


is calculated by means of the two methods presented above for different alignment 
geometry and triaxiality. For this task one chosen the most common quasiparticle 
configuration of one proton and one neutron with 7 = 7’ = 11/. The results of the 
exact diagonalization procedure and that obtained with the Schrédinger equation 
are visualized in Fig.3. First of all, one must remark the strong effect of even a small 
amount of tilting. Indeed, the degeneracy of the partner states achieved at high spins 
when the quasiparticle alignments are axial, is transformed into an increasing 
energy shift as a function of total angular momentum, regardless of the method. 
The exact calculations are overestimated in the chiral vibration region represented 
by the decrease of the energy difference with spin, and then underestimated in the 
higher spin regime, where few distinct dynamical scenarios can occur depending 
on the tilting. As the most experimental data are encompassed in the first vibrational 
regime, a suitable adjustment of the scale 1/% can in principle achieve an 
approximate equivalence between the two approaches. Comparing the whole 
spectrum, one observes that the minimum in the energy difference which marks the 
end of the chiral vibration regime is shifted to higher angular momentum values in 
comparison to the exact diagonalization results. The discrepancy between two 
calculations at higher spins increases for larger tilting angles. Nevertheless, the 
Schr6édinger equation results exhibit the same qualitative behavior as a function of 
total angular momentum as the exact problem. 
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(a) a=a’'=0° 
Schréd. eq. 
Schréd. eq. /-1 

— Exact diag. 


AE [1/4] 


(b) a=ar=1° 
Schréd. eq. 
Schréd. eq. J-1 

— Exact diag. 


AE [1/4] 


Schréd. eq. 
Schréd. eq. /-1 
— Exact diag. 


8 10 12 14 #16 18 #20. 22 
I 
Figure 3: The energy splitting of the chiral bands with j = 7’ = 11/2 as a function of 
angular momentum, determined within the Schrédinger equation and the exact 
diagonalization for y = 90°, 100° and tilting angles a = a’ = 0° (a), 1° (b), 3° (c). 
The Schrédinger equation results shifted to 7 — 1 position are also given as a 
reference. 
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A close inspection of Fig.3 suggests that the two calculations are congruent within 
a positive translation of angular momentum d/ ~ 1 in the case of the Schrédinger 
problem. More precisely, the results of the exact diagonalization for an angular 
momentum J is equivalent with the results for the Schrédinger equation 
corresponding to J+ 1. This relationship can be traced back to the dequantization 
and quantization procedures leading to the Schrédinger picture, where the 
uncertainty principle changes the original meaning of the angular momentum [46]. 
The correspondence J @ J + 1 between the two results is indeed one to one for the 
most part of the angular momentum range. Only the initial low spin part slightly 
deviates from this rule for more axially symmetric deformation (y 2 90°). As the 
diagonalization of the Schrédinger equation is significantly more computationally 
demanding than the exact diagonalization of the original cranked rotor 
Hamiltonian, their relationship can be used to determine model parameters of the 
first more easily. The usefulness of the more intricate Schrédinger picture resides 
in the separation of the kinetic and potential energy for a continuous variable 
associated with the total angular momentum projection, which allows the 
phenomenological interpretation of the quantum states in terms of anharmonic 
oscillations of the total angular momentum vector [38, 46]. The absolute energies 
visualized in Fig.4, are quite different in the two approaches. The excitation 
energies in the ground and excited bands obtained within the Schrédinger equation 
approach agree with those of the exact diagonalization mostly for very triaxial 
deformations y = 90° + 5° and low angular momentum states. The ground band 
levels calculated in the space of angular momentum states are underestimated up to 
T= 13 and /= 12 for axial and a = a’ = 3° quasiparticle alignments when maximal 
triaxiality y = 90° is assumed. Above these levels the exact excitation energies begin 
to be increasingly underestimated. This relationship depends on the triaxial 
deformation, such that at sufficiently axial deformation (y > 90°) all exact 
diagonalization excitation energies become underestimated. This is the result of the 
contradictory evolution of the energy spectrum from maximal triaxiality to more 
axial deformation. Indeed, Fig.4 (a) and (c) shows that the Schrédinger equation 
spectrum becomes more compressed with the increase of y, while the exact 
diagonalization spectrum become slightly expanded at low spins. In comparison to 
the ground band, the exact diagonalization results for the levels of the excited band 
are all underestimated. Excited energy levels of the exact diagonalization increase 
with the y deformation, whereas those coming from the Schrédinger equation have 
a similar behaviour only for axial alignments. When tilting is present, the excited 
states from the Schrédinger equation behave quit erratically at more axial 
deformations. As can be observed in Fig.4 (d), the Schrédinger equation spectrum 
for large values of y, is elevated in energy at low spins and lowered in energy at 
higher spins. A more general comparison between the axial and non axial alignment 
results, reveals the fact that the Schrédinger equation approximates better the 
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ground band for axial alignments and the excited band for tilted alignments. 


100 


100 


(a) s=1, a=a’=0° — Schréd. eq. + (b) s=2, a=a’=0° — Schréd. eq. 
tf gel Exact diag. { tog. Exact diag. - 4 


E,(1)-E1(8) [1/0] 
Eo(I)-E1(8) [1/40] 


100 ———— —— 
+ (c) s=1, a=a’=3° + —— Schréd.eq. j 
Foe Exact diag. 


E,(/)-E,(8) [1/0] 


—  Schréd. eq. 
----- Exact diag. 


Figure 4: Energy spectra of the ground band for axial (a) and tilted (c) alignments, 
and of the excited band for axial (b) and tilted (d) alignments, as a function of 
deformation y. The absolute energies are calculated for j = /’ = 11/2 and given in 
units of 1/% in reference to the ground band level J = 8. 


In order to demonstrate the practical usefulness of the Schrédinger equation 
approach, one considers here a particular numerical application performed on the 
observed partner bands in the !*°Pm nucleus [50]. Previous descriptions of the chiral 
bands in this nucleus performed with the Schrédinger equation approach [37, 38], 
considered maximal triaxiality y = 90° and established that one of the quasiparticles 


Schrédinger Equation for the Chiral Geometry in Atomic Nuclei 21 


is axially aligned, while the other has a sizable ~ 8° tilting. Note that due to the fact 
that at y = 90°, two MOI are equal (A; = Az), the problem with 7 = j’ cannot 
distinguish which quasiparticle is tilted, as the two alignment possibilities become 
equivalent. A microscopic investigation of the two chiral partner bands must be 
invoked in order to make an informed choice for the tilted quasiparticle. The 
findings of Ref.[51] can perfectly serve this role, which established that the 
dominant quasiparticle configuration in the structure of the chiral partner bands in 
'38Pm is represented by a proton particle hi1/2 spin with a substantial deviation from 
the s axis. Therefore the tilting in the present model for this nucleus is associated 
with the quasiparticle spin j’. The energy difference between the two partner bands 
obtained form fitting experimental data for '*Pm, is confronted with the 
corresponding experimental quantity in Fig.5. The model reproduces correctly the 
evolution with spin of AZ, including the decreasing part associated with the chiral 
vibration, the increasing part corresponding to the stabilized rotations, and the 
correct position of the minimum marking the transition between the two distinct 
regimes. Fig.5 also shows the energy difference resulting from the exact 
diagonalization of the original Hamiltonian containing angular momentum 
operators with the same scaling, alignment and deformation parameters. As was 
previously pointed out, these results shifts the minimum to lower angular 
momentum but qualitatively presents the same behaviour. It is worth to emphasize 
the fact that the evolution of the observed energy difference cannot be reproduced 
with considering tilting. 
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Figure 5: The evolution of the experimental energy difference between partner bands 
of the '*8Pm nucleus [50] as a function of total angular momentum is compared with 
the results of Schrédinger equation fits as well as with associated exact 
diagonalizetion results. The theoretical results correspond to J) = 35.68 MeV"! and 
tilting angles a = 0°, a’ = 8°. 
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Beside the energy levels, it is also instructive to compare the wave-functions 
resulting from the two formalisms. This is done in Fig.6, where one compared the 
coefficients of the total wave function expansion in rotation matrices obtained from 
exact diagonalization of the cranked rotor Hamiltonian and those resulted from the 
continuous density probability distribution offered by the Schrédinger equation. 
The results are shown only for the states with experimental chiral partners. It is 
remarkable to find that the two approaches are very much alike, with a predilection 
for higher K values in case of the exact diagonalization expansion coefficients. The 
ground band states have a single peak initially localized around K = 0 and which 
moves to positive K values with the increasing of total spin. On the other hand, the 
chiral partner states exhibit a node separating two peaks. This structure is well 
known to be associated with vibrational excitations, where the two peaks 
correspond to the turning points of the classical oscillation. The interpretation of 
this regime as vibrational cannot be confidently deduced from the discrete 
distribution of the K components, but it is naturally emerging from the continuous 
density distribution coming from the Schrédinger equation. This is one of the 
advantages of the Schrédinger equation over the solution of the problem in discrete 
angular momentum states. However, as spin increases, one of the two peaks starts 
to disappear by feeding the other. The surviving peak is opposite to the one from 
the ground band. One can therefore affirm that there is a transition from a chiral 
vibration to a static rotation where the two partner states achieve opposite 
handedness. 

The phenomenological analysis of the states resulting from the Schrédinger 
equation is completed by plotting in Fig.7, the chiral potential and the effective 
mass associated with angular momentum states considered in Fig.6. The transition 
between vibrational regime and the stable rotation mode discussed above can now 
be better understood by correlating it with the evolution of the chiral potential. The 
potential have a single asymmetric minimum at the lowest spin state and then starts 
to develop a second x < 0 minimum higher in energy. At higher spins, the second 
well completely traps the excited state forcing a stable rotation with K < 0 on it. In 
what concerns the symmetric mass function, it evolves from a flat shape towards a 
more parabolic x dependence as angular momentum is decreased. 


Schrédinger Equation for the Chiral Geometry in Atomic Nuclei 23 


Aix |? 


Aix |? 


Aix |? 


Aix |? 


[Aix |? 


K(x) K(x) 


Figure 6: The total wave function components Ajx of the first (s = 1) and second (s 
= 2) solutions spread over K = -J, -J + 1, ..., J- 1, J, determined from the exact 
diagonalization of equation (2.6) (blue histogram) and deduced from the 
Schrédinger equation (3.10) (red circles) when y = 90° and a= 0°, a'= 8°, matching 
the experimental data fits for the '*’Pm nucleus. The continuous line represents the 
density distribution p;(«) = |Fe(x)|? in respect to the dx metric. 
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Figure 7: The evolution with total angular momentum of the quantum chiral 
potential V(x) (a) and of the B(x) mass function (b), obtained from the numerical 
application on '38Pm nucleus for the states 7= 11, 12, 13, 14 and 15. The first two 
excited states, s = 1 (right) and s = 2 (left), are visualized for each angular 
momentum relative to their corresponding potentials. 


5. Conclusions 


The chiral geometry in atomic nuclei was modeled by means of a rigid alignment 
hypothesis applied to a Particle-Rotor Model extended to two quasiparticles. The 
rigid or frozen alignment approximation renders the relevant part of the Particle- 
Rotor Hamiltonian to a form of a simple rotor for the total angular momentum 
constrained by cranking terms resulting from the quasiparticle alignments. The 
resulted quantum problem is solved in the space of angular momentum states. 
Alternatively, one constructed a Schrédinger equation for a continuous projection 
variable, by means of a semiclassical description of the same cranked Rotor 
Hamiltonian. The two approaches are compared to each other in regard to the 
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corresponding energy spectra, in order to ascertain the deformation and alignment 
conditions for their equivalence. It is found, that although the results are 
quantitatively different, the spectral properties and different dynamical regimes of 
the chiral partner bands are qualitatively similar in both cases. A possible 
correspondence of the energy difference between the partner bands calculated in 
the two approaches is identified and proposed for the optimization of numerical 
applications. The two formalisms were also confronted in what concerns the 
properties of their corresponding wave function. This is realized by taking as a 
concrete example, the chiral partner bands observed in the '*8Pm nucleus. It thus 
found that both models provide the same evolution of the wave function with the 
total angular momentum, which is transiting from a vibrational regime towards a 
stable rotation mode where the partner bands have well established opposite 
handedness. In conclusion, the Schrédinger equation is a good approximation for 
the exact diagonalization results for low and intermediate spin values when the 
core’s shape is very triaxial. Otherwise, the Schrédinger equation approach can be 
considered as a stand alone model with similar spectral features and dynamic 
behavior to those of the exact problem. 
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